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ABSTRACT 

We present spectroscopic tests of MIR to FIR extinction laws in IRDC G028.36+00.07, a potential 
site of massive star and star cluster formation. Lim & Tan (2014) developed methods of FIR extinction 
mapping of this source using Spitzer-MlVS 24/rm and Herschel-PACS 70/im images, and by comparing 
to MIR Spitzer-IRAC 3-8/im extinction maps, found tentative evidence for grain growth in the highest 
mass surface density regions. Here we present results of spectroscopic infrared extinction (SIREX) 
mapping using Spitzer-IRS (14 to 38 ^m) data of the same IRDC. These methods allow us to first 
measure the SED of the diffuse Galactic ISM that is in the foreground of the IRDC. We then carry 
out our primary investigation of measuring the MIR to EIR opacity law and searching for potential 
variations as a function of mass surface density within the IRDC. We find relatively flat, featureless 
MIR-FIR opacity laws that lack the ^ 12 /im and 35 /im features associated with the thick water 
ice mantle models of Ossenkopf & Henning (1994). Their thin ice mantle models and the coagulating 
aggregate dust models of Ormel et al. (2011) are a generally better match to the observed opacity 
laws. We also find evidence for generally flatter MIR to FIR extinction laws as mass surface density 
increases, strengthening the evidence for grain and ice mantle growth in higher density regions. 
Subject headings: ISM: clouds — dust, extinction — infrared: ISM — stars: formation 


1. INTRODUCTION 

Massive star and star cluster formation are key pro¬ 
cesses for understanding galaxy evolution, since massive 
stars are vital sources of feedback for regulating the in¬ 
terstellar medium (ISM) and further star formation ac¬ 
tivity, while star clusters are the basic building blocks of 
galactic stellar populations. Despite this importance, we 
understand massive star and star cluster formation only 
very poorly because the regions forming such objects are 
relatively rare and thus typically far away (at least a few 
kpcs for most Galactic sources) and also deeply embed¬ 
ded inside dense molecular clouds. 

A better understanding of massive star and star clus¬ 
ter formation should result if we can find and charac¬ 
terize the initial conditions of these processes. Two key 
methods to achieve this involve studying the amount of 
IR extinction and sub-mm/mm emission of dust grains 
in cold, dense regions of molecular clouds in order to 
measure mass surface densities and thus masses of cloud 
structures (see, e.g. Tan et al. 2014, for a review, here¬ 
after TI4). To do this accurately, one needs to know 
the opacities and emissivities of the dust grains, includ¬ 
ing potential systematic evolution in these properties. 
This is the goal of our study. A better understanding of 
dust properties will also be important for modeling heat¬ 
ing and cooling in the clouds, astrochemistry, and the 
ionization fraction, which affects dynamics by enabling 
coupling of the mostly neutral gas to magnetic fields. 


The NIR to FIR properties of dust extinction (from I 
to 25 //m) have been measured in the Lupus molecular 
cloud using background stars by Boogert et al. (2013). 
They found evidence for water ice (presumably forming 
on dust grains) on sightlines with Ay > 2 mag and that 
with increasing Ak the > 5 pm continuum extinction in¬ 
creases relative to Ak- The most extincted line of sight 
that they probed had Ak = 2.46 mag {Ay 20 mag). 
These are relatively low mass surface density conditions 
compared to the clouds that are thought to form massive 
stars and star clusters (T14). Lutz et al. (1996) stud¬ 
ied the MIR extinction law (from 2.5 to 9 pm) towards 
the Galactic Genter using hydrogen recombination lines. 
They found a relatively flat extinction law, e.g., com¬ 
pared to bare grain models of Draine (1989). This sight¬ 
line has Ay ^ 30 mag, but is essentially probing diffuse 
ISM conditions that are not likely to be close to initial 
conditions for star formation. McGlure (2009) studied 5 
to 20 pm dust extinction properties on sightlines through 
molecular clouds with 0.3 < Ak < 7 towards background 
stars, finding systematic evolution with increasing Ak 
that was attributed to ice-mantle mediated grain coag¬ 
ulation (see also Flaherty et al. 2007; Ghapman et al. 
2009; c.f., Ascenso et al. 2013, who did not find evidence 
for extinction law variations up to Ay 50 mag). 

Infrared dark clouds (IRDCs) (e.g., Perault et al. 1996; 
Garey et al. 1998; Ormel et al. 2005; Rathborne 
et al. 2005; T14) cast shadows at MIR wavelengths, 
10 pm. They are occasionally opaque at far-infrared 
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wavelengths, up to ~ 100 ^m, against the Galactic dif¬ 
fuse background emission (e.g., Lim & Tan 2014, here¬ 
after LT14) due to their very high mass surface densities 
~ 1 g cm“^ {Ay ~ 200 mag) and cold < 15 K temper¬ 
atures. Being massive clouds with such physical condi¬ 
tions, IRDCs are likely to be representative of the initial 
conditions for massive star and star cluster formation 
(T14). 

Butler and Tan (2009, 2012; hereafter BT09, BT12) 
studied mass surface densities, E, of 42 starless cores in¬ 
side 10 IRDCs via MIR extinction (MIREX) mapping 
by utilizing Spitzer-IRAC 8 /im data of the GLIMPSE 
survey (Ghurchwell et al. 2010). Butler, Tan and Kain- 
ulainen (2014, hereafter BTK14) used deeper Spitzer- 
IRAC 8 /im images toward one of the 10 BT09 IRDGs 
(IRDC G028.37-b00.07, also referred to as IRDC C), in¬ 
cluding NIR extinction offset corrections of Kainulainen 
and Tan (2013), to produce a higher dynamic range ex¬ 
tinction map. 

The advantages of such MIREX mapping are high an¬ 
gular resolution (2"), a near-fully sampled image (except 
in regions of bright MIR emission) and an independence 
to dust temperature (except, again, in very warm regions 
that are sources of MIR emission). A major uncertainty 
is the need to estimate the background intensity by inter¬ 
polation from regions around the cloud. The background 
intensity can also have significant small scale fluctua¬ 
tions, which however can be assessed from the off-cloud 
regions, but which are an inherent random uncertainty. 
Another major issue is the need to estimate the level of 
the MIR foreground intensity from the diffuse ISM, since 
IRDCs are typically at several kpc distances. The key 
method to estimate the foreground is to measure it em¬ 
pirically in spatially independent “saturated” regions of 
the IRDC where essentially negligible background light 
makes it through the cloud (as long as such regions exist), 
and then assume it is spatially invariant across the IRDC 
or some local region (BT12). The accuracy of measuring 
the saturated intensity level depends on the photometric 
precision of the image and sets a maximum Egat that can 
be probed in the map. 

MIREX mapping does require knowing the ^ 8 /rm 
opacity per unit mass: different dust models, such 
as the thin and thick Ossenkopf and Henning (1994, 
hereafter OH94) moderately coagulated models or the 
Weingartner & Draine (2001) bare grain models, show 
only modest variation at the level of ~ 30%. Then 
a (refractory)-dust to gas mass ratio also needs to be 
assumed (1:142) is our adopted fiducial value (Draine 
2011; note BT09 and BT12 assumed 1:156). Over¬ 
all, BT09, BT12 and BTK14 adopted the OH94 thin 
ice mantle model (10^ years of coagulation at density 
^ 10® cm“®) as a fiducial, which has (averaged 

over the IRAC band) ~7.5 cm^ g“^. With these meth¬ 
ods, BT12 could trace details of highly extincted re¬ 
gions up to, Esat ~ 0.5 g cm“^, i.e., Ay ~100 mag. 
BTK14 could reach somewhat higher Egat ~ 1 g cm“^, 
i.e.. Ay ~ 200 mag. 

LT14 extended MIREX methods to the FIR regime 
by using Spitzer-MIPS 24 pm MIPSGAL data (Carey 
et al. 2009) and Herschel-PACS 70 pm archival data 
(proposal IDs: KPGT-okrause-1 & KPOT-smolinar-l) 
of IRDC G028.37-1-00.07. From studying ~ 1' regions 



Fig. 1.— Upper: IRDC G028.37+00.07 (Cloud C) as viewed 
by Spitzer-YRAG 8 fim (blue), Spitzer-MIPS 24 pm (green) and 
Herschel-PACS 70pm (red). The dashed ellipse shows an approxi¬ 
mate boundary for the IRDC defined by Simon et al. (2006). The 
strip observed by Spitzer-IKS is shown: Order 1 (19.5—38.0pm) was 
measured in the white rectangle region; Order 2 (14.0—21.3 pm) in 
the blue rectangle. Order 1 and 2 data are both available where 
these regions overlap, i.e., in the magenta/violet region, which has 
been divided into 21 10^^ x 10^^ regions (numbered from upper 
right), from which individual spectra have been analyzed. Magenta 
squares indicate “dark” regions that are used for our extinction 
law study. Violet squares are regions we consider contaminated by 
24 pm point source emission. The two green regions at the outer 
ends of the strip are the regions used to assess the spectrum of the 
background (BG) emission from the diffuse ISM. Small white circles 
indicate the positions of the dark, saturated cores, Cl, C4, Cll, 
that were investigated by LT14. Bottom: Mass surface density 
map (in g cm“^) derived from Spitzer-MIPS 24 pm data, 
of the same region as the RGB map above. The slit positions of 
Orders 1 and 2 are shown. 

around three dense, locally saturated cores (the white cir¬ 
cles in Figure 1), LT14 found tentative evidence of k{X) 
variation with E in the sense that the relative extinction 
curves appear to flatten as E increases. Such behavior 
is consistent with an evolution from the OH94 thin ice 
mantle model to the thick ice mantle model, in which 
there is a growth of the volume of the ice cover from 
0.5 X that of the refractory component to 4.5x. Some 
degree of flattening could also result from coagulation of 
the grains (Cardelli et al. 1989; Weingartner & Draine 
2001; Ormel et al. 2011, hereafter Oil). 

In this paper, we present a new spectroscopic IR ex¬ 
tinction (SIREX) mapping method by using Spitzer-IRS 
mid-infrared spectroscopic data that covers a part of 
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IRDC G028.37+00.07. This IRDC (which we will also 
refer to as Cloud C, BT09) is located at 5 kpc dis¬ 
tance and is one of the most massive IRDCs known with 
~ 7 X lO"* Mq within an effective radius of about 8 pc 
(BTK14) (i.e., the dashed ellipse shown in Figured]). It 
hosts a modest level of star formation (e.g, the MIR- 
bright sources in Fig. dl see also Zhang et al. 2009), 
but overall the star formation activity appears to be rel¬ 
atively low, and so this cloud should represent an earlier 
stage of massive star cluster formation. 

The main objective and plan of the paper are as 
follows. The SIREX method adopts saturation based 
MIREX (BT12) and FIREX (LT14) methods to then al¬ 
low investigation of the extinction law variation as a func¬ 
tion of mass surface density. In §2 we review and update 
the results of LT14 for extinction law variation in dif¬ 
ferent mass surface density regions, including utilization 
of WISE 12 p,m data. We introduce the detailed meth¬ 
ods of SIREX mapping in §3. In §4, we present our re¬ 
sults, first demonstrating the measurement of the Galac¬ 
tic background and foreground from the diffuse ISM, and 
then presenting the extinction law in different mass sur¬ 
face density environments. We discuss the implications 
of these results and concluded in §5. 

2. PHOTOMETRIC ANALYSIS 

Here we update the photometric MIREX and FIREX 
analysis methods of LT14 as applied to IRDG C, in¬ 
cluding: utilizing the improved 8 /rm extinction map of 
BTK14, which probes to higher values of E; comparing 
to the DL07 diffuse ISM dust emission model; analyzing 
WISE 12 /xm data; and comparing to Oil dust opacity 
models, in addition to those of OH94. 

2.1. Dust Model Opacities in MIR to FIR Photometric 

Bands 

Following the methods of LT14, we consider several 
different dust models from OH94 and now also Oil. 
The thin & thick ice mantle models of OH94 (ice 
composition of H 20 :CH 30 H:GO:NH 3 = 100 : 10 : 1 : 1 ) 
and the models of Oil (ice composition of 
H 20 :CH 30 H:CO:NH 3 = 100 : 66 : 5 : 5 ) are moderately 
coagulated grains (i.e., coagulation for various periods 
at riH = 10® cm“® for OH94 and at nn = 10® cm“® for 
Oil. For consistency with LT14, we utilized 10® yr of 
coagulation OH94 models (which is equivalent to 10® yr 
of coagulation at nu = 10® cm“®). We also adopt the 
gas to (refractory component) dust mass ratio of 142 
(Draine 2011) that is used in LT14. 

Densities of nn ~ 10® cm“® are typical of dense cores 
and clumps within IRDCs (e.g.. Tan et al. 2014), how¬ 
ever the lifetimes of these structures are less well con¬ 
strained. We expect the clouds to be at least one local 
free-fall time old, i.e., 1.4 x 10® yr at this density, but 
ages that are ^ 10 x longer are certainly a possibility 
and may help account for high levels of deuteration that 
are seen in some IRDCs (Kong et al. 2015). The longest 
duration Oil model involves coagulation for 10^ yr. The 
applicability of such timescales to dense clumps and cores 
seems less likely to be valid, especially given the appar¬ 
ent turbulent nature of these molecular clouds, which 
should induce density fluctuations on shorter timescales. 
We note that the OH94 and Oil are local models, valid 


for particular, uniform densities (c.f., the structured core 
models of Weidenschilling & Ruzmaikina 1994). We note 
also that the OH94 and Oil models involve dust grain 
coagulation that is mediated by relative grain velocities 
that are set by a model of subsonic turbulence. On the 
larger clump scales within IRDCs, we expect supersonic 
turbulent motions to be more relevant, although these 
may be moderated by the presence of strong magnetic 
fields (Filial et al. 2015). Supersonic turbulence is seen 
to decay within a few free-fall times in numerical simula¬ 
tions (Stone et al. 1998; Mac Low et al. 1998; see review 
by McKee & Ostriker 2007), so on the smaller scales of 
dense cores we may expect conditions to have evolved 
closer to a subsonic turbulent cascade. 

We derive filter and background spectrum weighted 
opacity values. BT09 derived these for IRAC bands 1- 
4 and the MIPS 24 /xm band using the diffuse Galactic 
background spectrum from Li & Draine (2001; hereafter 
LDOl) and a total to refractory dust component mass 
ratio of 156. LT14 modified these values based on the 
total to refractory component dust mass ratio of 142 and 
expanded to Herschel-PACS 70, 100 and 160 /xm filters. 
Here, we recalculate these opacities by utilizing the dif¬ 
fuse Galactic background emission model of DL07, and 
now also including the WISE 12/xm band filter (Tabled]). 

These opacities can be used to calculate E from a given 
set of imaging data. Their relative values will also be 
compared to those observed in IRDG G. 

2.2. The MIR to FIR Dust Extinction Law in IRDCs 
from Broad Band Photometry 

We now re-visit the tentative evidence presented by 
LT14 of extinction law variation, i.e., flatter k \/in 
higher E regions of IRDG C. We utilize the same MIR 
and FIR images processed by LT14, i.e., Spitzer-IBAC 
3.5, 4.5, 5.9, 8/xm, Spitzer-MlPS 24^m and archival data 
of Herschel-PACS 70 /xm. The angular resolutions of the 
IRAC bands, MIPS 24 /xm band and PACS 70 /xm band 
are ~2", 6" and 6", respectively. 

We also calibrate and use WISE 12 /xm imaging 
data with angular resolution ~6.5". Following MIREX 
(BT09, BT12) and FIREX mapping (LT14) methods, we 
estimated the background emission by applying the small 
median filter (SMF) method, i.e., a median filter size set 
to 1/3 of the major axis of the IRDC, i.e., 4', in the re¬ 
gion outside the dashed ellipse shown in Fig. |T] followed 
by interpolation inside the ellipse. 

We also adopted the saturation based foreground esti¬ 
mation method of BT12: the values of the foreground, 
relative to that in IRAC band 4 at 8 /xm, were mea¬ 
sured towards cores Cl, C4 and Cll and are shown in 
the top panel of Fig [2] They are compared to the DL07 
model of Galactic plane diffuse ISM emission (excluding 
stars), which is shown by the solid line, with the diamond 
shaped boxes showing convolution of this model with the 
instrumental filter response functions. 

The Eg^m map of IRDC C from BTK14, which probes 
to higher mass surface densities than the BT12 map, is 
now utilized so that we are able to consider binned val¬ 
ues up to 0.4 g cm“^, compared with the 0.3 g cm“^ of 
LT14. Following LT14, we derive optical depth maps in 
1' regions around each saturated core Cl, C4 and Cll in 
each waveband and then maps of relative opacity, nor¬ 
malized to that of IRAC band 4 at 8/xm, and with angular 
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Fig. 2.— (a) Top: Total IR to sub-mm SED of the Galactic plane diffuse ISM, normalized to the Spitzer-IKAC 8 pm band flux. The 
solid line shows the dust-only component of the DL07 model and the dotted line indicates the total flux from dust and stellar components. 
Large diamonds indicate convolution of the dust-only SED with the filter response functions of corresponding instruments (panel b). The 
red, blue and green open circles show foreground intensities, relative to the IRAC 8 pm value, measured towards the three saturated cores 
Cl, C4 and Cll, respectively, (b) Second from top: Filter response functions of IRAC bands 1-4 (solid lines), WISE W3 band at 12 pm 
(dotted), MIPS 24 pm (dashed) and PACS 70, 100 Sz 160 pm (dot-dashed), (c) Third from top: Opacities relative to effective opacity 
measured in the IRAC 8 pm band. Black dotted/dashed/dot-dashed lines show thin/thick ice mantle and bare grain dust models of OH94 
after 10^ yr of coagulation at = 10® cm“^ (or 10® yr of coagulation at = 10® cm“®). Blue, green, gold and red solid lines indicate 
Oil fully ice coated dust models after 0.1, 1, 3 and 10 Myr of coagulation at = 10® cm“^, respectively. The red/green/blue filled 
squares (for C1/C4/C11 regions; slight wavelength offsets applied for clarity; see also zoom-in insets) show K,u/K,8iJ,m for the indicated 
ranges. For uncertainties in these data, see Figure [3] (d) Bottom: Same as (c), but now showing opacities relative to effective opacity 
measured in the MIPS 24 ^m band. 
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TABLE 1 

Telescope Band and Background-Weighted Dust Opacities Per Total Mass^ (cm^ g~^) 


Dust Model*^ 

IRAC3.5 
3.52 pm^ 

IRAC4.5 
4.49 fj.m 

IRAC6 
5.91 pm 

IRAC8 
7.80 /im 

WISE12 
12.0 pm 

MIPS24 
23.0 /j,m 

PACS70 
74.0 fim 

PACSIOO 
103.6 pm 

PACS160 
161.6 pm 

OH94 thin mantle 

26.9 

(19.5) 

14.4 

(11.7) 

10.31 

(9.37) 

8.19 

(8.04) 

12.8 

6.26 

1.14 

0.598 

0.287 

UH94 thick mantle 

52.4 

(38.0) 

17.5 

(14.3) 

20.1 

(18.2) 

10.3 

(10.1) 

36.7 

11.4 

3.96 

1.26 

0.399 

on 30ky 

12.3 

(11.9) 

7.03 

(6.91) 

4.79 

(4.75) 

7.14 

(7.12) 

8.24 

5.01 

0.79 

0.38 

0.16 

on O.lMy 

16.3 

(14.1) 

8.16 

(7.57) 

5.14 

(4.97) 

7.60 

(7.56) 

9.87 

6.59 

1.22 

0.590 

0.242 

on 0.3My 

41.4 

(25.7) 

18.4 

(12.2) 

9.08 

(6.64) 

8.91 

(8.13) 

11.5 (11.0) 

7.34 (7.29) 

1.25 

0.59 

0.24 

on IMy 

115 

(61.3) 

70.5 

(37.8) 

40.3 

(22.0) 

21.7 

(14.1) 

22.5 (16.0) 

11.4 (9.17) 

1.34 (1.19) 

0.57 (0.53) 

0.21 (0.20) 

on 3My 

30.9 

(17.9) 

31.5 

(17.7) 

32.4 

(17.8) 

26.9 

(15.7) 

25.4 (15.7) 

20.1 (12.8) 

7.30 (4.36) 

3.49 (2.09) 

1.13 (0.70) 

on lOMy 

10.8 

(6.65) 

7.22 

(4.55) 

5.53 

(3.57) 

5.05 

(3.59) 

5.41 (4.02) 

4.74 (3.49) 

4.45 (2.89) 

4.84 (2.93) 

5.14 (2.87) 


^ A total to refractory dust mass ratio of 142 is adopted (Draine 2011). OH94 opacities in the IRAC bands have been scaled from values in 
parentheses to include an approximate estimate of the contribution from scattering, based on the Weingartner & Draine (2001) dust model 
for Ry = 5.5 dust. The contribution of scattering to the total extinction is expected to be very small for A > 10 /i.m for these dust models. 
The online Oil models provide scattering opacities, so we compute both the pure absorption extinction (in parentheses; shown when there is 
a significant difference) and total extinction including scattering. 

^ OH94 models for 10^ yr of coagulation at uh = 10^ cm~^ (equivalent to 10® yr of coagulation at nn = 10® cm”®). Fully ice covered Oil 
models for various coagulation times as noted (at nn = 10® cm”®). 

^ Mean wavelengths weighted by filter response and background spectrum. 



0.1 0.2 0.3 0.4 0.5 

(9 ) 




Fig. 3. Scatter plots of ^24/im/^8/Aini ^70^im/^8/^m5 ^8/im/^24/im; Kl2/jm/^24/^m and K70/jm/^24/^m verSUS 

measured from the BTK14 Ss/im map and imaging data from Spitzer-IHAC 8 fim, WISE 12 /im, Spitzer-MlPS 24 pm and Herschel-PACS 
70 pm. The size of pixels in all maps have been regridded to 6" resolution (i.e., the approximate beam sizes of the 12, 24, and 70 pm 
images). Red, green and blue dots correspond to pixels in Cl, C4 and Cll regions, respectively. The open squares are the mean values of 
logio of relative k, with error bars indicating the Icj dispersion. The black open squares are the overall average of the three cores. Black 
dotted, dashed and dot-dashed lines are relative k. values of OH94 thin, thick ice mantle and bare grain models, respectively. Blue, green, 
gold and red solid lines indicate Oil fully ice coated models of 0.1, 1, 3 and 10 Myr of coagulation time, respectively. 
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resolution of ^ 6", set by the PACS 70 //m image. We 
also consider opacities normalized to those derived from 
MIPS images at 24 /rm. 

Note, that this analysis, following BT12 and LT14, as¬ 
sumes that the IRDC itself is a negligible source of emis¬ 
sion, either intrinsically (e.g., from thermal emission of 
surface layers of warm dust, including transiently heated 
small grains) or from scattering of radiation into our line 
of sight from other directions by dust within the IRDC. 
This assumption appears to be supported by the fact 
that independent dark saturated regions are seen within 
the IRDC across the studied MIR to FIR wavebands. If 
scattering has an effect, it is expected to be most im¬ 
portant at the shorter, MIR wavelengths, where it has 
been seen to occur in deep IRAC 3.5 mum and 4.5 /rm 
images of nearby cores, producing the “core shine” phe¬ 
nomenon (e.g., Pagani et al. 2010; Lefevre et al. 2014). 
Accurate assessment of the effects of such scattering in 
IRDCs requires radiative transfer modeling and assump¬ 
tions about the radiation fields and 3D IRDC structure. 
Given our focus in ths paper on dust extinction prop¬ 
erties at longer wavelengths, especially extending from 
^ 10 /rm to ^ 70 /rm, where scattering should be of rela¬ 
tively minor importance for most grain models, we con¬ 
tinue with the approximation adopted in our previous 
studies of neglecting the effects scattering. 

Each pixel in the maps has a value of Es/im- Scat¬ 
ter plots of ^ 12 /im/ ^24:fini! ^8/im and ^70^in/^8/im 
versus are shown in the top row of Fig. [3l The 

bottom row shows the equivalent plots of Ksfim/K 24 fi,m, 
Ki2tJ.m/K24fj.m and Kro^tm/ii24tj.m versus Eg^m- Note, with 
the new BTK14 E map, including its NIR offset correc¬ 
tions, there are relatively few pixels with E < O.lg cm“^. 
Thus, we replace the low-E bin of 0.1 ± 0.05 g cm“^ used 
by LT14, with a bin defined by E = 0.15 ± 0.05 g cm“^. 
Binned averages for Eg^m = 0.15 ± 0.05,0.2 ± 0.05, 0.3 ± 
0.05, 0.4 ± 0.05 g cm“^ are shown - these are values of E 
that should not be affected by saturation in the BTK14 
map (these values are also shown in the third and fourth 
panels of Fig. [2|). Higher values of binned Eg^ni are 
shown in Fig. [H but these results should be treated with 
caution since they may be affected by saturation, espe¬ 
cially affecting both 8 and 24 /rm results. 

Comparing to the models of OH94, the observed 
«^ 70 /im/Kg^m and K 70 iim/K 24 tj.m values are closer to the 
thick ice mantle case than the thin ice mantle case (sim¬ 
ilar to the findings of LT14). However, the 12 /im WISE 
results do not favor the thick ice mantle model (see also 
Fig. [21 Comparing to the Oil models, the observed 
Krofim/ and urofim/K 24 fj,ni values are consistent with 
models with coagulation times of >3 Myr. However, 
again the shorter wavelength results (including the IRAC 
bands 1 to 3; see Fig. El) tend to favor the shorter coagu¬ 
lation time models. In other words, there is not a single 
dust evolution model that provides a fully consistent fit 
to the data. 

We next examine whether there is any trend of evolving 
opacity law with E. Inspecting the panels of Fig. El we 
notice that the /C 70 Mm/K 8 /jm and /C 70 Mm/K 24 Mm values ap¬ 
pear to increase as E increases. We examine the Pearson 
correlation coefficients, r, and p-values (probability of 
chance correlation) for these cases finding r = 0.25, 0.34 
and p = 5 X 10“®, 1 x 10“® when analysing all the data 


for K 70 Mm/Kg^m and /t 70 ^m/K 24 /xm, respectively. These 
results appear to indicate evidence for a systematically 
evolving MIR to FIR extinction law, i.e., becoming flat¬ 
ter, as mass surface densities (and presumably also vol¬ 
ume densities) increase. Such a flattening is expected 
from models of grain growth (e.g., the Oil models). The 
other wavelength data do not show trends with E that 
are as significant. However, this may be expected since, 
e.g., the Oil models also do not predict as significant a 
variation as they do in comparison with the evolution of 

We conclude that the 70 /im data in comparison with 
the 8 and 24 /xm data show evidence for growth of dust 
grains via coagulation (Oil) and/or thickening of ice 
mantles (OH94, although such models struggle to match 
shorter wavelength results) in the densest regions. How¬ 
ever, the particular shape of the extinction law and its 
possible variation are only quite poorly constrained. 

3. SPECTROSCOPIC IR EXTINCTION (SIREX) MAPPING 
METHODS 

3.1. MIR Spectroscopic data of IRDC G028.37+00.07 

We analyze archival Spitzcr-YRS Long-Low spectra 
(program: 3121, PI: Kraemer, K. E.) with spectral res¬ 
olution R ~57-126 of IRDC G028.37-P00.07 (Cloud C). 
The IRS Long-Low module consists of two long slits (1st 
Order and 2nd Order; hereafter LLl and LL2) that have 
a size of 168"xl0.5" each. LLl and LL2 observe the 
wavelength ranges of 19.5-38.0 pm and 14-21.3 pm, re¬ 
spectively. The IRS spectral data of IRDC C include 4 
Long-Low observations of staring modes that cover ~ T 
in length for each order (LLl and LL2) diagonally across 
the IRDC (Fig. |T]). The overlapped region of LLl and 
LL2 is ^ 4' long. The extraction of the observed spectra 
is executed by utilizing the CUbe Builder for IRS Spec¬ 
tra Maps (CUBISM, Smith et al. 2007). The extracted 
spectra before any corrections are shown in Figure E] 

3.2. Radiative Transfer Assumptions 

We follow the same simple one-dimensional radiative 
transfer model of the MIREX (BT12) and FIREX (LT14) 
methods and here review the basic equations that de¬ 
scribe the input quantities that are required for the anal¬ 
ysis. We need to know the intensity of radiation directed 
toward the observer at the location just behind, /i,p, and 
just in front, 1^,1, of the target IRDC. The infrared emis¬ 
sion from the IRDC (both intrinsic emission and that 
scattered into the line of sight) is assumed to be negligi¬ 
ble so that 

Iu,i = e ^"1^,0, (1) 

where optical depth is total (absorption + 

scattering) opacity at frequency n per unit total mass and 
E is the total mass surface density. The spectrum of 
i.e., the background, is to be estimated from calibrated 
spectra of “off-positions” away from the target IRDC 
region, while I^^i is derived from the observed spectrum 
toward the location of interest. 

However, we cannot observe I^p and I^^i directly be¬ 
cause these intensities are contaminated by foreground 
emission, R,foT:e so that the observed background emis¬ 
sion, Iup,ohs, and the observed intensity towards the 
cloud, /y,i,obs, are 

obs — Ii>P F 7^/fore 
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Fig. 4. — Extracted spectra of the selected 14 “dark” regions 
(i.e., without contamination of point sources) arranged in order of 
increasing S 24 ^m (shown in top left corner of each panel in g cm~^; 
the number in parentheses is the region number, see Fig. 1). The 
summed spectrum from the background (BG) regions is shown in 
the lower-right panel. The photometric error bars are estimated 
by CUBISM. All spectra are shown before the zodiacal correction 
has been applied. 

and 

obs (^) 

respectively. Here we are assuming that Iu,iore, which is 
emitted from the diffuse ISM, does not vary significantly 
in its intensity over the face of the IRDC and to the 
off-position. 

Therefore, as with MIREX and FIREX mapping, 
SIREX mapping also requires knowing the intensity of 
the foreground emission toward the target region. While 
the MIREX and FIREX derived Iu,iore from observed 
“saturated” regions, i.e., they have same observed in¬ 
tensity that agrees with the minimum intensity within a 
range set by 2a (or other set range) noise level of the im¬ 
ages, SIREX mapping toward IRDC C cannot apply this 
same saturation-based foreground method, since there is 
no slit position that is coincident with a saturated core, 
i.e., as defined from the MIPS 24 /rm maps (LT14) (see 
Fig. 1). Thus we need to apply another method to esti¬ 
mate the foreground spectrum. 

3.3. Background and Foreground Estimation 

We select two background regions, i.e., off-positions, 
from near the north (BGl) and south (BG2) ends of the 
slit position (^50"xl0" on each side). Unfortunately, 
these regions are not covered by both orders (Fig. [T]). 
The selected background regions have a ~I2% difference 
in their average Spitzer-MlPS 24 fim photometric inten¬ 
sities, <5% differences in Spitzer-IRAC 8/xm, ~I0% dif¬ 
ferences WISE 12 fj,m intensities and ~35% differences 
Herschel-PACS 70 pm. These gradients are utilized, be¬ 
low, for linear interpolation of the background spectrum 



Fig. 5.— (a) Top: Spitzer-IKS spectrum from the background 
(BG) regions (green rectangles in Fig. la). The black solid line 
shows the originally extracted spectrum without any correction 
(hereafter BG raw) with photometric errors. The black dotted line 
is a model spectrum of ~0.8x the fiducial zodiacal light model 
(hereafter zodi). The blue solid line is the background spectrum 
after being corrected, i.e., BG raw - 0.785x zodi. (b) Middle: BG 
region spectra showing the effect of different fractions of zodiacal 
light subtraction as a correction to the background spectrum. The 
constraint of matching the Spitzer-M.IPS 2A{im filter-weighted data 
point, i.e., the average of the small median filter background model 
of LT14 of the two BG regions (red diamond) and also of 21 regions, 
sets our choice of using 0.785 x the fiducial zodi model intensity. 
The green and blue diamonds indicate similar Spitzer-IKAC 8 pm 
and WISE 12 pm photometric data, respectively, for the BG re¬ 
gions. (c) Bottom: Spectral line masked (see text), smoothed (to 
~ 0.5 pm resolution at ~ 24 pm) and zodi-corrected background 
spectrum (dashed line), and then iteratively de-reddened via OH94 
thin ice mantle model (see text) background spectrum (solid line). 
If the 011 1 Myr dust model is used for de-reddening the back¬ 
ground then the result is the dotted line, i.e., showing about 10%— 
20% differences depending on wavelength. The dot-dashed line 
indicates the DL07 model of the diffuse Galactic background emis¬ 
sion. 

along the slit. 

We compare the Spitzer-IRS spectroscopic data with 
the Spitzer-MIFS 24 fim image of the background and 
selected 21 regions. We subtract fiducial model zodiacal 
emission from the extracted spectra to obtain intrinsic 
intensities as a function of wavelength. However, we find 
that the intensities of the MIPS 24 fim filter-weighted 
IRS background spectrum after these zodiacal light cor¬ 
rections have ^10% higher values than in the MIPSGAL 
image. In order to correct this offset, we consider dif¬ 
ferent fractions of the fiducial model zodiacal light sub¬ 
traction from the spectra of the background and the 21 
selected regions (Fig. [5l middle panel). We find that the 
average value of the best matched fraction of the zodiacal 
model from all regions is 78.5%. 

While the BGl and BG2 regions are near the ends 
of the IRS long-slit aperture, they are still quite close 
to the IRDG and are thus likely to be afi'ected by 


















Lim, Carey & Tan 




Fig. 6. — SED of the diffuse Galactic ISM. (a) Top: Derived 
SED of the background emission, 1^,0, from the BG regions. The 
red solid line shows the de-reddened and foreground-subtracted 
Spitzer-IKS data (see text). The small diamonds show photometric 
estimates from the BG regions (based on SMF background models). 
The black solid line shows the DL07 model (dust only; best-fit 
normalized to the Spitzer-IKS data) and the large diamonds show 
its convolution with the filter response functions of the Spitzer- 
IRAC bands, WISE 12 /rm, Spitzer-y\TPS 24 /rm and Herschel- 
PACS 70 pm bands, (b) Bottom: As for (a) but now showing 
derived SED of the foreground emission towards IRDC C. Now the 
small diamonds show photometric measurements towards the three 
saturated cores Cl, C4 and Cll. 

extinction from the IRDC and its surrounding CMC. 
Thus we de-redden the BG spectra to obtain the in¬ 
trinsic background spectrum that is not contaminated 
by the IRDC. We first flag and mask the molecular 
emission lines (S(III):18.7I3 /rm, H 2 S(0): 28.221 ^m, 
S(III): 33.480 ^m and Si(II): 34.815 /im) to simplify the 
spectrum. Then, de-reddening is applied to estimate 
via ^ djy fore)o ^ (hig. [HI)- 

To estimate the amount of de-reddening, we adopt the 
^ 24 tin\ map from the work of LT14 (see Fig. 1) and k^, 
from the OH94 thin ice mantle model. Note, the va¬ 
lidity of this assumption can be checked after the fact 
by comparing our derived BG spectrum with the DL07 
(empirically-calibrated) model, but, still, this is a fun¬ 
damental limitation of the SIREX method when the off 
positions are too close to the IRDC, which ultimately 
introduces a systematic uncertainty in the shape of the 
derived IRDC extinction curves, although not affecting 
their relative shapes. To test the sensitivity of our results 
to this choice, we also try de-reddening the background 
emission using the Oil model of 1 Myr coagulation time. 
This leads to a ~ 20% higher intensity of the derived 
background (at 24 fim) than when using the OH94 thin 
ice mantle model (Figure [5]). 

Having estimated the background spectrum, we uti¬ 
lize it to estimate the foreground spectrum by scaling 
the background by a factor, ~0.35, derived by consid¬ 
ering the ratio of the average MIPSGAL 24 /im inten¬ 
sity towards the three saturated cores Gl, C4 and Gil 
compared with the average MIPSGAL 24/im SMF back¬ 
ground map intensity of the BGl and BG2 regions. 

This estimated foreground is then used, iteratively, in 
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Fig. 7.— Derived background spectrum (dashed line; grey 
shaded region shows estimated errors — see text), foreground spec¬ 
trum (scaled up by a factor ~0.35 from the background), and 
zodi-corrected observed spectra (solid lines) from each of the 14 
24 gm-dark regions (separate panels, in order of increasing S 24 ;im 
(indicated in g cm“^ in top left corners, along with region numbers 
in parentheses). 

the above de-reddening analysis, since the foreground 
is not extincted by the IRDC. We find this iterative 
method converges quickly, yielding our final estimated 
background spectrum (/i/,o) from the BG regions and 
scaled foreground spectrum (/j/,fore)- These are shown in 
Figure [HI together with our photometric estimates that 
probe these spectra. We also compare to the best-fitting 
DL07 model. 

After obtaining an estimate of the background spec¬ 
trum from the BG regions, we then consider how it may 
vary spatially along the IRS strip. To do this we uti¬ 
lize WISE 12 /im, Spitzer-M\PS 24 /tm and Herschel- 
PACS 70 /im photometric estimates of the backgronnd 
intensities (via the SMF method) at each end of the slit. 
Spatial gradients at these wavelengths are measured: at 
12 and 24 /im the intensities vary by about 10% across 
the length of the slit; at 70 /im the variation is about 
35%. We use these results to carry out linear interpola¬ 
tion of the background intensity as a function of position 
and wavelength. The estimated background, i.e., /,y,o.obs, 
is then compared to the observed emission, i.e., I^^i^ohs, 
from each of the 21 dark regions along the slit length. We 
find that an overall scaling of a 15% increase in the back¬ 
ground intensity is necessary to achieve self-consistent 
opacity results (see S]). This correction is within the 
estimated size of systematic errors associated with the 
above methods, and may be due to some combination of 
errors in the background de-reddening, background spa¬ 
tial interpolation, foreground fraction estimation and zo¬ 
diacal correction methods. Onr final adopted estimates 
for /jy^o.obs and I^^i^ohs are shown for each region in Fig¬ 
ure [T] 





















9 



Fig. 8. — Extracted spectra from the 10”xl0” dark regions (col¬ 
ored lines) and the foreground/background model (dotted/solid 
black lines). The comparison between the foreground spectrum 
and spectra of S 24 ^m ^0.25 g cm~^ regions indicates we see satu¬ 
ration at A <17/rm. Triangles indicate 24/rm photometric data of 
the saturated cores Cl, C4, Cll that help to normalize the fore¬ 
ground model (also shown in|6]l. Diamonds are Spitzer-MlPS 24/rm 
photometric data points for all regions, which helped normalize the 
zodiacal subtraction. 

The uncertainties of are constrained based on the 
observed color variations over the 12/im and 24/ini maps, 
the observed MIPS 24 /im intensity fluctuations in the 
surrounding regions outside the IRDC C ellipse, and the 
intrinsic photometric errors (that were shown in Fig.|4]). 
By considering the scatter in pixel by pixel values of 
h 2 /im/l 24 /im in an elliptical annulus 0.7 to 1.3x the semi¬ 
major axis of the IRDC, we derive a la uncertainty, i.e., 
dispersion, of about 7% in the spectral slope of our esti¬ 
mate of the background spectrum. The color variation of 
the wavelength range between 24 and 38 /im is uncertain 
and so is potentially larger. Thus, to be conservative, 
we assume an uncertainty in the spectral slope in this 
wavelength range of of 20%. The fluctuations in 
l 24 tin\ are also estimated in the same annulus and found 
to have a Icr dispersion of about 18%, which we adopt 
for our error analysis. Thus the overall uncertainties in 
Ij/.o, including photometric errors, are summed and are 
shown in the panels of Fig.[7]as the grey shaded regions. 
These panels also show the scaled foreground and the 
zodi-corrected observed spectra for each 24 /im-dark re¬ 
gion. These spectra are also shown together in Figure [H 
These are the required quantities to estimate the MIR to 
FIR extinction law. 

4. RESULTS 

4.1. Spectrum of the Galactic plane diffuse ISM 

Figure [6] shows that our estimates of the spectrum of 
the SED of the diffuse ISM generally agree well will the 
DL07 model. Our derived spectral slope of this SED 
from the Spitzer-IRS data rises somewhat more steeply 
towards longer wavelengths than the DL07 model. This 
could either indicate a real difference or may be due to 
systematic errors introduced in our analysis, such as the 



X (iim) 


Fig. 9.— Upper panel: The optical depth (t>) curves of the 
14 selected dark positions of the Spitzer-IRS slit on the IRDC 
G028.37-1-00.07. The color of spectra are same as in Fig. |8]and are 
listed in the legend of the lower panel. Lower panel: Relative k 
curves (normalized to K 24 ^im) of the same regions are compared to 
model relative opacities of OH94 thin (black dotted line) and thick 
ice mantle models (black dashed line), OH94 bare grain model (dot- 
dashed line) and Oil models of 0.1, 1, 3 and 10 Myr of coagulation 
time (blue, green, gold and red solid lines). 

de-reddening method, which assumed the OH94 thin ice 
mantle model opacities. 

As an independent check on our derived spectrum of 
the diffuse ISM, we can search for evidence of saturation 
at the short wavelength end of the observed spectra in 
high E regions. Figure Hand [5] shows evidence for such 
saturation in the four highest E regions (all with Eg/im > 
0.2 g cm“^) at A < 18 /tm. Note the normalization of the 
foreground model (dotted lines in this figure) is based on 
the 24/rm photometry towards the three saturated cores 
Cl, C4 and Cll, so if our model of emission from the 
diffuse ISM were in error, then we would not expect to 
see saturation, i.e., coincidence of the foreground model 
with the observed spectra. 

4.2. The MIR to FIR Extinction Law and Evidence for 
Grain Growth 

We calculate the optical depth, of each selected 
10"xl0" 24 /rm-dark region via Ti,=ln(/i,_o/7i/,i), and 
these are shown in the upper panel of Fig. H As ex¬ 
pected, higher E 24 ;_im regions show higher optical depths. 
Using E 24 ;im (from the LT14 FIREX map), we then de¬ 
rive Ki, = E24;im/'r;y, and show K\lK24i_im in tbe lower 
panel of Fig. H along with the OH94 and Oil dust mod¬ 
els. The absolute values of k\ are shown in Fig. [TUI 

To quantify the shape of the extinction laws that are 
probed by the Spitzer-IRS data we fit various spectral 
indices, /3, where KxfK 24 fim oc A^). We define Index A 
as Pa being measured from 23 to 30 pui; Index B as 
Pb being measured from 30 to 34 /rm; and Index C as 
Pc being measured from 23 to 34 /rm. The divisions of 
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Fig. 10.— The opacity (fty) curves (thick solid grey line; shaded areas show uncertainties) of the individual regions (each panel lists 
S 24 |jm in g cm“^ and region number in parentheses), compared to model kx of OH94 thin ice mantle (dotted), thick ice mantle (dashed) 
and bare grain (MRN) model (dot-dashed). The Oil models of 0.1, 1, 3 and 10 Myr of coagulation time (blue, green, gold and red solid 
lines) are also shown. 


A and B at 30 /im were chosen to accentuate the “V”- 
shaped water ice feature that is prominent in the thick ice 
mantle models of OH94. Figure [IT] shows these various 
indices and the combination of ,5^ — 13a as a function 
of S 24 ^m- The uncertainties of the spectral indices of 
each region are approximately estimated by adopting the 
range of indices that are allowed from the upper and 
lower limits of the SIREX-derived opacities (i.e., from 
the shaded areas of Figure ITUl). 

The absolute values of these indices and the combi¬ 
nation f3B — Pa are consistent with the OH94 thin ice 
mantle model. The values of /3 b disfavor the thick ice 
mantle model. Oil models with ^ 1 Myr of coagulation 
do reasonably well, but in general the Oil models are 
too flat. However, we caution that there are systematic 
errors (e.g., the dereddening of the background regions) 
that could be affecting the derived spectral indices of the 
extinction laws. 

Next we look for systematic variation of the extinc¬ 
tion laws with mass surface density. All three indices 
and the combination /3 b — Pa show hints of flattening 


as ^ 24 tin\ increases. The Pearson p values are reason¬ 
ably significant for Index B and C (~0.008 and ^0.01, 
respectively). Thus Index B shows relatively significant 
changes with increasing S 24 ^nij which may indicate that 
water ice is in fact beginning to change the shape of the 
extinction curves in the denser regions (although not to 
the extent of OH94 thick ice mantle models). Such an 
interpretation would help to reconcile the results of the 
Spitzer-IRS extinction curves with the photometric re¬ 
sults of lj2l which we present together in a summary in 
Figure [12 

5. CONCLUSIONS 

We have presented a new method of MIR to FIR 
spectroscopic extinction (SIREX) mapping of an IRDC, 
which can observationally trace grain growth processes 
in a region that is expected to form massive stars and 
a massive star cluster. This is the first study to deter¬ 
mine the detailed extinction law in the wavelength range 
14 pm to 38 pm in such dense regions. In combination 
with our MIREX and FIREX maps, we are also able to 
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Fig. 11.— Spectral indices of extinction laws versus S 24 ^m. (o.) 
Top left: “Index A,” the spectral slope, 13a (where KA/^ 24 f 3 m oc 
A^^), measured from 23 to 30 pm. The black circles show the 
results of the fit to the 14 analyzed 24 pm-dark regions, with best- 
fit relation shown by thick solid line. /3 a values of the OH94 and 
Oil dust models are marked with horizontal lines, (b) Top right: 
Same as (a), but now for “Index B,” which is evaluated from 30 to 
34 pm. Note the relatively large value of the OH94 thick ice mantle 
model, due to the water ice feature, (c) Bottom left: Same as (a), 
but now for “Index C,” which is evaluated from 23 to 34 pm. (d) 
Bottom right: fds — Pa versus S 24 ^m ■ 

search for systematic variation in the extinction law as a 
function of mass surface density. 

Our new results, especially from WISE 12 /tm MIREX 
mapping and Spitzer-]RS SIREX mapping, indicate that 


KA/N 24 /rm is relatively flat from 5 to 25 /im. The 
strong water ice absorption feature of the OH94 thick ice 
mantle models is not seen. SIREX results then indicate 
a relatively steep drop of kx/K 2 Aprm out to 30 ^m, but 
then hints of some flattening, especially in denser regions, 
out to 35 ^m. The FIREX results at 70 /im suggest 
this flattening extends out to these long wavelengths and 
is also stronger in higher density environments. Com¬ 
parison with OH94 and Oil dust models suggests that 
this may require some combination of water ice mantle 
growth and grain coagulation. 

One implication of these results is that there could 
be very significant grain growth on large scales within 
IRDCs, so that pre-stellar cores should be modeled with 
such large grains. This can potentially affect cooling 
rates (generally leading to cooler equilibrium tempera¬ 
tures; e.g., Keto & Caselli 2008). It can also affect astro- 
chemistry and ionization fractions. 

Another important implication is for mass estimates of 
dense regions from sub-mm and mm dust emission. The 
difference in derived mass surface density that results 
from using either thin or thick OH94 dust opacities is 
about a factor of three. 

Future studies of spectroscopic observation toward 
dense cores and clumps by using finer angular resolu¬ 
tion and higher sensitivity (e.g., JMdST-MIRI) will help 
to better explore the evidence of grain growth in different 
star-forming environments. 
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trophysics Data Analysis Program grants ADAP10-0110 
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Fig. 12.— Overall 24 /rm-normalized MIR to FIR extinction law, combining photometric (MIREX and FIREX) results of [[2] (from 
Fig.®) and spectroscopic (SIREX) results. The solid lines with black-white gradient indicate SIREX KA/^ 24 ^m results at different E 24 ^m- 





















